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Abstract. Black hole (BH) X-ray binaries (XRBs) are X-ray luminous binary systems comprising a BH accreting matter from 
a companion star. Understanding their origins sheds light on the still not well understood physics of BH formation. M33 X-7 
hosts one of the most massive stellar-mass BH among all XRBs known to date, a 15.65 M :, BH orbiting a 70Mr;i companion 
star in a 3.45 day orbit. The high masses of the two components and the tight orbit relative to the large H-rich stellar component 
challenge our understanding of the typically invoked BH-XRBs formation channels. The measured underluminosity of the 
optical component further complicates the picture. A solution to the evolutionary history of this system that can account for 
all its observed properties has yet to be presented, and here we propose the first scenario that is consistent with the complete 
set of current observational constraints. In our model, M33 X-7 started its life hosting a 85-99 Mm primary and a 28-32 M . 
companion in a Keplerian orbit of 2.8-3.1days. In order to form a BH of 15.65Mq, the initially most massive component 
transferred part of its envelope to the companion star and lost the rest in a strong stellar wind. During this dynamically stable 
mass transfer phase the companion accreted matter, to become the presently underluminous 70 star. 
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1. INTRODUCTION 

X-ray binaries (XRBs) are key systems for studying the physics of black holes (BH). In fact, although solitary BHs 
do not emit electromagnetic radiation, they become detectable X-ray sources when they have a stellar companion 
transferring matter to them. M33 X-7 provides a unique physical laboratory for the study of BHs, massive stars and 
XRBs. This system harbors the heaviest star ever discovered in an XRB orbiting one of the most massive BHs found in 
this class of systems, a 70M© star orbiting a 15.65 M© BH every 3.45 days. The massive components and the tight orbit 
mark M33 X-7 as an evolutionary challenge. The search for a plausible evolutionary scenario is further complicated 
by the luminosity of the stellar component, which is lower than what evolutionary models predict for a single star of 
70M©[1]. Our analysis aims to find a solution to the evolutionary history of M33 X-7 that can simultaneously explain 
all its observed properties. 



2. M33 X-7'S CURRENT PROPERTIES AND DISTANCE UNCERTAINTY 

In Tablel we summarize some of the observed system's properties for a distance to M33 as adopted by the discovery 
team, and we discuss how some of these parameters would change if a different distance is considered. In fact, M33's 
distance is not well constrained by the literature, and various studies place it between 750-1017kpc [2, 3, 1, 4, 5, 6, 
7, 8, 9, 10, 11, 12, 13, 14]. At the present time, the BH is accreting part of the matter the companion star is losing via 
stellar wind. In fact, given the extreme mass ratio, a phase of mass transfer (MT) through Roche-lobe (RL) overflow 
would be dynamically unstable and rapidly evolve into a merger of the two components. 



TABLE 1. Observed parameters for M33 X-7. The BH mass (Mbh), its compan- 
ion mass (M2), Spectral Type, effective temperature (T e ff), and luminosity (L2), the 
orbital eccentricity (e), and inclination (;'), and M33 distance (d), are taken from 
Orosz et al. [1]. Pietsch et al. [15] determined the orbital period (P). The X-ray 
luminosity (L x ) accounts for variations in the X-Ray flux over different observa- 
tions [16, 17, 15, 18, 1, 19]. If the full distance range of 750-1017kpc is adopted, 
via the ELC [20] we calculate M 2 ~ 55-103 M , M BH ~13.5-20M G , i ~77°-71°, 
log(L 2 /X ) ~5.62-5.89, andL x between -lxlO 37 and ~3xl0 38 erg sec" 1 . 



Parameter 


Value 


Parameter 


Value 


M BH (M q ) 


15.65 ± 1.45 


P (days) 


3.45301 ± 0.00002 


M 2 (Mq) 


70.0 ± 6.9 


e 


0.0185 ± 0.0077 


Spectral Type 


07 III to 08 III 


iO 


74.6 ± 1.0 


?eff(K) 


35000 ± 1000 


L x (10 38 ergs" 1 ) 


0.13 to 2.49 


Log(L 2 /L ) 


5.72 ± 0.07 


d (kpc) 


840 ± 20 



3. MODELING THE SYSTEM 

To investigate a solution to M33 X-7's history, we perform binary evolution calculations considering a variety of 
progenitor masses and initial orbital periods. Guided by the luminosity of the stellar component, lower than what is 
expected from evolutionary models of single stars, we explore binary sequences where the BH companion accretes 
mass from the BH progenitor. Furthermore, given the short orbital period, we consider binary systems that start their 
life already in a tight orbit, hence undergoing a phase of MT during the Main Sequence (MS) phase of the primary 
component. The stellar evolution models are calculated with an up-to-date version of Eggleton's stellar evolution code, 
STARS[21, 22, 23, 24]. Given the spatial metallicity gradient of M33 [2], we assume a metallicity 50% of the solar 
value for all our models. 



3.1. Scan of the Parameter Space 

We evolve binary systems from the Zero-Age MS until the end of the primary's MS considering progenitor masses 
between 20-1 30 M and 10-100M Q for the primary and the secondary, respectively, and initial orbital periods ranging 
from 1 to 10 days. We perform a first scan of the parameter space rejecting the sequences where the primary overfills 
its RL only after the end of its MS. Furthermore, given the present high mass of the star, we also exclude the sequences 
where the secondary transfers mass back to the primary after having accreted from it. In our model, assuming quasi- 
conservative MT (with the only mass loss from the system being due to the stellar wind of both components), the 
primary transfers most of its H-rich envelope to the secondary, and becomes a Wolf-Rayet star (WR). The stronger 
WR wind then interrupts the MT and blows away the remaining primary's envelope to expose its massive He core. 

At the end of the primary's MS, we use evolutionary models of single He stars (see Appendix A) to determine 
the mass lost from the system, and the consequent change in orbital period during the primary's core He burning 
phase, until collapse l . We then reject the sequences where, at He exhaustion, the components masses are lower 
than the minimum observed values. Since no other episodes of RL overflow can have occurred from the primary's 
collapse, we evolve each secondary as a single star, and reject the sequences where the mass of the star does not fall 
within the observed range when the model matches the currently observed effective temperature and luminosity. After 
the primary's collapse, we study the evolution in time of the orbital separation, eccentricity, and spin of the stellar 
component accounting for torque between the binary components (following Belczynski et al. [26], and references 
therein 2 ), changes in the stellar radius during the MS lifetime of the companion star, stellar wind mass loss [28], orbital 



The final stages of the primary's life during and beyond carbon burning last too short to significantly change the system's parameters (~60yr for 
an initially ~25M He star)[25] 

2 For the 2 nd order tidal coefficient £2 we use stellar models from Claret [27] for masses of ~ 63M.0 and ~ 79.5 Mq, and derive E2 = 
-5.4566 - 7.37243 ■ *ms° 562 > wnere 'MS is time in units of the MS lifetime. The formula has only a very weak dependence on the initial mass 
for the range of interest. 



angular momentum loss due to gravitational radiation [29], and accretion from the companion's stellar wind onto the 
BH [30]. We consider the BH as a point mass. We explore a variety of initial orbital configurations, by scanning the 
parameter space made up of the kick magnitude (Vt), orbital separation (a pos tBH), an d eccentricity (e pos tBH), following 
Kalogera [31], Willems et al. [32]. Specifically, we consider kicks between 0-1300 km/s, orbital separation between 
0-100R©, and eccentricities between and 1. We assume an isotropic distribution for the direction of the kick. We 
interrupt the calculation when the orbital period crosses the observed value and the eccentricity of the orbit falls within 
the observed range. Finally, of the sequences that fulfill all the above requirements, we reject the sequences where the 
star's radius at present is bigger than the distance from the center of the star to the LI point [33]. 



3.2. Example of a Successful Evolutionary Sequence 

In Figure 1 we show the results for one of the successful evolutionary sequences, typical of M33 X-7's history. 
The progenitors comprise a primary and a secondary of ~97M© and ~32M©, respectively, in an orbit of ~2.9days. 
During the first ~1.8 Myr the evolution is driven by mass loss via stellar winds, causing an expansion of the orbit to 
~3. 25 days. While still on its MS the primary overfills its RL and begins MT onto the secondary. This stronger mode 
of mass loss brings the primary out of thermal equilibrium, and in response the star shrinks, recovering its thermal 
equilibrium while always maintaining hydrostatic equilibrium, and hence keeping the MT dynamical stable. During 
the MT phase, while the more massive primary is transferring mass to the less massive secondary, the orbit shrinks, 
but when the secondary becomes the more massive component the orbit begins expanding [34]. The primary transfers 
most of its H-rich envelope and becomes a WR star, and the strong WR wind (<~2 to 3T0~ 5 M yr _1 ) eventually 
interrupts the MT. During the ~99,000 years of conservative MT, the primary becomes a ~5 1 M WR star, while 
the original 32 M secondary becomes a massive ~69M© O-type star. Once the WR wind sets in and the primary 
detaches, the wind blows away the remaining primary's envelope to expose the ~25M© He core. At the same time, 
the now more massive secondary is losing mass via its own wind at a lower rate (~ lO _6 M yr _1 ). This mass loss 
causes the orbital period to further increase until the end of the primary's MS and throughout its core He burning 
phase. At this time the orbit is circular and the spin period of each star is expected to be synchronized with the orbital 
period. At primary's collapse, after ~3.7Myr, M33 X-7 hosts a ~16M© evolved WR star, and a ~64.5M© O-star 
companion in a ~3.5day orbit. Unable to support itself through further nuclear fusion, the WR star collapses into a 
BH, and 10% of the rest mass energy is released as the BH's gravitational energy. Additionally, collapse asymmetries 
and associated neutrino emission may impart a kick to the newly born BH, even with no baryonic mass ejection at 
collapse. Both these effects modify the orbital configuration, slightly shrinking the orbit to ~3.4days, and inducing 
an eccentricity. For the remaining ~0.2Myr, the evolution of the system is driven by mass loss via the secondary's 
stellar wind, causing a further expansion of the orbit, and bringing the orbital period to the currently observed value. 
The fraction of this stellar wind accreted by the BH is too small to significantly influence the orbital evolution, but it 
is adequate to explain the observed X-Ray luminosity. At the present time, after ~3.9Myr, M33 X-7 comprises a BH 
of ~14.4M© and an underluminous O-star of ~64M© in a slightly eccentric ~3.45day orbit. 



3.3. Parameter Space Left 

All the successful sequences follow a path qualitatively very similar to the specific example described in detail 
here. Adopting a M33 distance of 840±20kpc, the sequences that match all observed properties within la errors, 
are constrained to host 96-99 M primaries, and 32 M (within 1M uncertainty) secondaries in orbits with initial 
periods of 2.8-2.9days. If the full distance range of 750-1017kpc is adopted, the progenitors are constrained to host 
primaries between 85-99 M , secondaries between 28-32M , and initial orbital periods between 2.8-3.1 days. Our 
model excludes secondaries at present more massive than ~65 M Q because they fail to explain the observed luminosity 
(see Figure 2). Two factors contribute to the apparent underluminosity of the stellar component (see appendix B). On 
one hand, the orientation of the system with respect to our line of sight and associated projection effects reduce 
the star's measured luminosity (~87%); on the other hand the secondary was not born as a ~63-65M© star, but 
instead accreted much of its mass from the BH progenitor (~13%). Looking at the system's properties at primary's 
collapse (Figure 3), our models show no ejection of baryonic mass. Furthermore, the allowed eccentricities post-BH 
formation are small (0.012-0.026). However, given the lack of kinematic information, we can not exclude kicks as 
high as ^850km/s. This apparent discrepancy is explained by the change in orbital inclination at BH formation, which 




FIGURE 1. Evolution of the stellar and orbital parameters of M33 X-7. From the top: components masses (M), secondary's 
effective temperature [log(r e ff)] and luminosity [log(L2)], orbital period (P), and eccentricity (e). From the left, the different 
evolutionary stages are: the beginning of the MS (a), the MT phase (b), the end of the MS (c), the core He burning phase for 
the primary (d), and the post BH formation phase (e) until the present time (note the non uniform x-axis). The grey-shaded areas 
represent the observational constraints as reported in Table 1. For the case shown, a kick of 120 km/s is imparted to the newly born 
BH. 




FIGURE 2. Current luminosity. The circles and triangles are the results of detailed binary star evolution calculations for all 
successful sequences for a M33 distance of 840 ± 20kpc, and of 750-1017kpc, respectively. The BH X-Ray luminosity (left), and 
the secondary's luminosity (right) are shown as a function of the secondary's mass at present. The dark and light grey shaded areas 
represent the observational constraints for a distance of 840 ± 20kpc, and of 750-1017kpc, respectively. Lx is calculated according 
to Bondi and Hoyle [30], and the error bars account for the uncertainties in the stellar wind parameters, and depict the highest and 
lowest Lx values. Some of the data points are omitted for clarity. 



increases with the kick magnitude. Hence, our model constrains the kick to point mostly orthogonal to the orbital 
plane, without significantly affecting the eccentricity. 




FIGURE 3. Orbital parameters at BH formation. The data points are as in Figure 2. (a) eccentricity as a function of the orbital 
separation post-BH formation; (b) mass of the BH progenitor and (c) change in the orbital inclination at BH formation as a function 
of the kick magnitude. The BH progenitor mass accounts for the 10% of rest mass energy released as the BH's gravitational energy 
at collapse. Given that the BH progenitor mass lies within the observed range for the BH mass, no baryonic mass is ejected at 
collapse. Some of the data points are omitted for clarity. 



4. CONCLUSIONS 

For the first time, we present a quantitative theoretical model that succeeds in explaining all observed properties of 
M33 X-7. The massive BH and its massive, but underluminous companion in a tight, slightly eccentric orbit are all 
physically explained, given our current understanding of massive binary evolution. Our results further shed light on the 
effect of mass transfer on massive stars and their internal structure, the relation between BH and progenitor masses, 
and asymmetries associated with BH formation through massive-star core collapse. 
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APPENDIX A 

Using the code STARS we create 3-25 M© He star models with the same input physics used for the binary evolution 
calculations. We evolve each model after the end of the core He burning phase and calculate the duration (fn e ) of 
this phase and the amount of mass lost (AM) in stellar wind as a function of the initial He star mass (M He ,i)- From 
these models we derive f He = 0.323221 +6.24256 -M^ 5762 where t He is in Myr, AM = 0.282 - 0.28 -M H e,i + 0.052 • 

M^ ;- 0.00102 -M^ j if M He ,i < 17M®, and AM = -23.5355 + 10.262 ■]n(Ar H e,i) if M He i > 17 M , where M He ,i and 
AM are in M©. 



APPENDIX B 

Correction to the luminosity due to the inclination of the system: The shape of the star in M33 X-7 is distorted by 
rotation and tides, which cause the equatorial regions to be colder than the poles (Von Zeipel theorem), and dimmer. 
Given that our stellar models do not incorporate these effects on the star's surface temperature (T e s) and luminosity (L), 
we use the ELC code to determine an appropriate correction to apply to our models, and we find that T e g — > 7^/0.954, 
and log(L) — > log(L) — 0.14. The observed luminosity is lower than the true luminosity because we are looking at the 
colder equatorial regions of the star. 

Correction to the luminosity due to the partial-rejuvenation of the secondary: Previous studies on the effect of 
accretion onto massive MS stars have shown that accreting stars not always rejuvenate, and that the result of mass 
accretion might be a star with a chemical structure unlike that of a single star of the same mass [35]. To a large extent, 
this effect is controlled by the assumed semiconvective mixing efficiency which, in turn, depends on the criterion 
for convection used in the stellar model. Braun and Langer [35] used the Ledoux criterion and showed that the non- 
rejuvenated models appear to be underluminous for their new mass during the remaining Main-Sequence evolution. 
Following this idea, we use the Ledoux criterion in the single and binary evolution calculations, after having calibrated 
the semiconvective efficiency parameter against some of the results reported by Braun and Langer [35]. Our results 
confirm that the rejuvenation of the secondary star after mass transfer is at most partial. 



